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SUMMARY 


Accurate  experimental  data  are  given  on  the  shape  and  the  location 
of  detached  shock  waves  on  cones  and  spheres  at  Mach  numbers  from  1.17 
to  1.81.  The  data  are  correlated  to  obtain  equations  that  describe  the 
shock  waves.  This  knowledge  of  the  shock  waves  should  be  useful  in 
calculations  of  the  pressure  distribution  and  the  pressure  drag  of  the 
fore  part  of  cones  and  spheres.  The  experimental  data  on  shock  waves 
are  compared  with  theory. 


INTRODUCTION 


One  of  the  problems  in  aerodynamics  that  is  of  current  interest  is 
the  problem  of  predicting  the  shape  and  the  location  of  detached  shock 
waves,  or  bow  waves.  The  problem  is  of  practical  importance  because  a 
knowledge  of  the  shock  wave  is  useful  in  calculating  the  drag  of  the 
nose  of  the  body  behind  the  shock  wave.  The  flow  behind  a  detached  shock 
wave  is  subsonic  in  a  region  that  is  bounded  by  the  shock  wave,  the  sonic 
line,  and  the  fore  part  of  the  body.  If  these  boundaries  of  the  subsonic 
region  are  known,  then  the  pressure  distribution  over  the  fore  part,  or 
nose,  of  the  body  and  the  pressure  drag  of  the  fore  part  of  the  body  can 
be  calculated  (references  1  and  2) . 

A  few  theoretical  attacks  on  the  problem  of  describing  the  detached 
shock  waves  on  bodies  of  revolution  have  been  made  (references  3  to  6). 
Attempts  have  been  made  to  compare  the  theoretical  results  with  the 
available  experimental  results  (references  3,  5,  and  6).  Experimental 
results  have,  however,  been  scarce. 

The  purpose  of  the  present  paper  is  to  give  experimental  data  on 
the  location  of  detached  shock  waves  on  spheres  and  conical-tip 
cylinders.  (For  brevity,  the  conical-tip  cylinders  are  referred  to 
in  the  present  paper  as  cones.)  For  this  purpose,  series  of 
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interferograms  were  taken  of  the  flow  around  spheres  and  cones,  of 
various  vertex  semiangles  between  35°  and  90°,  at  Mach  numbers  ranging 
from  1.17  to  1.8l.  Measurements  of  the  shock  waves  were  made  on  these 
interferograms  and  were  correlated  to  obtain  equations  that  describe 
the  part  of  the  shock  wave  ahead  of  the  sonic  line . 


c0->  cl*  c2>  •  •  • 


C 

d 

K 

M 


n 


P 

R 


S1 


s2 


S 

T 


SYMBOLS 

coefficients 

constant  of  apparatus 

diameter  of  sphere;  diameter  of  cylinder 
curvature  of  bow  wave  at  its  nose 

free -stream  Mach  number 

number  of  points  measured  on  a  bow  wave 

pressure  of  air  in  jet 

gas  constant 

axial  distance  of  shock  wave  from  center  of 
sphere 

axial  distance  of  shock  wave  from  plane  of. 
shoulder  of  cone 

fringe  shift  in  free  stream 

static  temperature  of  air  in  jet 

stagnation 'temperature  of  air  in  jet 
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coordinates 

dimensionless  coordinate  (x/d) 
dimensionless  coordinate  (y/d) 
ratio  of  specific  heats 
root-mean-square  deviation 
density  of  atmospheric  air 
density  of  air  in  jet 
vertex  semiangle  of  conical  tip 


APPARATUS 


A  3-inch-square  free  jet  and  an  interferometer  with  4-inch-square 
plates  were  used.  This  apparatus  has  been  described  in  reference  7- 
Different  Mach  numbers  were  obtained  by  using  five  supersonic  nozzles 
similar  to  those  described  in  reference  7-  For. the  calibration  of  each 
nozzle,  interferograms  were  taken  that  showed  part  of  the  free  stream 
and  also  part  of  the  undisturbed  region  outside  the  jet.  The  fringe 
shift  S  in  the  free  stream  was  determined  by  comparing  such  an  inter - 
ferogram  with  a  no-flow  interferogram.  Then  the  free-stream  density  p 
was  computed  by  means  of  equation  (3)  of  reference  7> 


p 

—  =  CS  +  1 
pa 


where  C  is  a  constant  of  the  apparatus  and  where  pa,  the  density  of 
the  undisturbed  air,  was  determined  from  the  measured  atmospheric 
pressure  and  temperature.  The  supply  to  the  jet  was  regulated  in  such 
a  way  that  the  jet  emerged  with  its  static  pressure  equal  to  that  of 
the  atmosphere.  *  Then,  the  static  temperature  T  in  the  free  stream 
was  computed  by  the  general  gas  law, 
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The  stagnation  temperature  Ts  was  measured  with  a  thermocouple  a  short 
distance  upstream  of  the  nozzle  where  the  airspeed  was  low.  Then  the 
Mach  number  was  obtained  from  the  relation 


Ts 

—  =*  1 
T 


+ 


7-zJmS 

2 


The  calibration  yielded  the  values  1.17,  1*30,  1-37,  1.62,  and  1.81 
for  the  Mach  numbers  of  the  five  nozzles.  The  determination  of  the 
Mach  number  given  by  each  nozzle  was  made  as  carefully  as  was  practicable. 
Some  uncertainty,  nevertheless,  exists  in  the  final  result.  This 
uncertainty  is  estimated  to  be  not  greater  than  3  percent. 


RESULTS 

Spheres 


Interferograms  of  the  flow  at  several  Mach  numbers  around  spheres 
of  various  diameters  d  are  shown  in  figures  1  to  5-  Additional 
interferograms  with  spheres  of  other  diameters  were  taken  at  each  Mach 
number.  At  a  given  Mach  number,  corresponding  distances  on  different 
interferograms  were  found  to  be  proportional  to  d;  that  is,  in  these 
tests,  the  Reynolds  number  had  no  apparent  effect  on  the  part  of  the 
shock  wave  that  was  ahead  of  the  sonic  line.  Therefore,  only  one 
interf erogram  is  shown  for  each  Mach  number  except  for  Mach  number 
1.81,  for  which  interferograms  of  the  flow  around  spheres  of  three 
different  diameters  have  been  included  to  illustrate  this  result  (fig.  5) • 

Axial  distance  of  shock  waves.-  The  axial  distance  sp  of  the  shock 
wave  from  the  center  of  the  sphere  has  been  obtained  from  figures  1  to  5* 
The  variation  of  the  quantity  sp/d  with  Mach  number  is  shown  in 
figure  6.  The  equation  of  the  curve  that  has  been  drawn  through  the 
experimental  points  is 


2  =  1  (M  -  D-1/3  k  (i) 

d  3 

Over  the  range  of  Mach  number  used,  the  curve  fits  the  data  well. 
Equation  (l)  is  not  applicable,  however,  at  Mach  numbers  much  greater 
than  the  largest  used  in  the  present  tests,  because,  according  to 
equation  (l),  the  distance  between  the  nose  of  the  shock  wave  and  the 
center  of  the  sphere  is  less  than  the  radius  of  the  sphere  at  Mach 
numbers  greater  than  3*37* 
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Whether  equation  (l)  remains  applicable  at  Mach  numbers  less  than 
the  smallest  used  in  the  present  investigation  is  not  yet  known.  The 
equation  does  give,  however,  the  correct  asymptotic  value,  infinity, 
for  the  quantity  sq/d  as  the  Mach  number  approaches  unity. 

The  theory  of  reference  3  predicts  the  variation  of  the  distance 
between  the  nose  of  the  shock  wave  and  the  nose  of  the  sphere, 

S1  1 

,  with  Mach  number,  for  Mach  numbers  near  unity.  A  comparison 

d  2 

of  the  theoretical  results  with  the  present  experimental  results  is 
given  in  figure  7*  One  of  the  curves  in  figure  7  is  taken  from 
figure  2  of  reference  3*  The  curve  for  the  experimental  results  is 
plotted  from  equation  (l)  of  the  present  paper.  Reference  6  also 
predicts  the  location  of  the  shock  wave.  One  of  the  curves  in  figure  7 
is  based  on  the  results  of  reference  6. 

Equations  of  shock  waves. -  Equations  for  that  part  of  the  shock 
wave  which  bounds  the  subsonic  region  have  also  been  determined.  The 
coordinates  x  and  y  of  points  on  the  shock  wave  were  measured  with 


series 

X  =  C0  +  c^Y  +  C2Y2  +  c3Y3  +  C^Y11  +  c5Y5  +  cgY6  +  .  .  . 

to  describe  it.  Because  of  the  choice  of  the  origin  and  because  the 
shock  wave  is  symmetrical  about  the  x-axis. 


c0  -  C1  c3  _  c5 


.  .  -  0 


t 
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and,  therefore, 

X  =  c2Y 2  +  c^  +  CgY6  +  .  .  . 


(2) 


The  first  two  terms  on  the  right-hand  side  of  the  equation  are  found  to 
he  sufficient  for  an  accurate  description  of  the  shock  waves. 

Application  of  the  method  of  least  squares  yields  the  following 
formulas  for  the  coefficients: 


■<=2  -  T  (XVl2!1!8  -  X*V  X1!6) 

=■*  -  xd^xv  -  x*v  xv) 


where  i  =  1,2,3,  .  .  •  n,  the  values  of  and  Y^  are  taken  from 

table  I,  and 


Associated  with  each  pair  of  values  C2  and  c^  is  the  quantity 


oh.  — 

where  5i  =  c2Yf +  c^Y^  -  X^.  The  value  of  5  is  a  measure  of  how  well 
the  equation  fits  the  data;  the  smaller  5  is,  the  better  the  fit. 

Computed  values  of  c2,  c^,  and  5  are  listed  in  table  II.  In 
figure  8,  these  values  of  c2  and  c^  are  plotted  against  M. 
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An  inspection  of  the  figure  shows  that  there  is  a  fairly  regular  variation 
of  and  with  Mach  number.  Curves  that  show  the  approximate 

variation  have  been  arbitrarily  drawn  in  figure  8  to  aid  in  interpolation 
at  other  than  test  Mach  numbers.  The  equations  of  these  curves  are 


c2 


1  _  e-0.70(M-l) 


(3) 


C4  =  -0.042(M  -  1) 


W 


The  scatter  in  the  values  of  c^  is  much  larger,  on  a  relative  basis, 
than  the  scatter  in  the  values  of  C2*  This  large  uncertainty  in  the 
value  of  c)|  does  not  necessarily  indicate  a  large  uncertainty,  how¬ 
ever,  in  the  values  of  X  calculated  from  equations  (2),  (3),  and  (4), 
because  the  term  in  equation  (2)  that  contains  c^  is  usually  much 
smaller  than  the  term  that  contains  C2- 


Curvature  of  shock  waves.-  Reference  5  is  a  theoretical  investigation 
of  the  problem  of  the  detached  shock  wave  in  front  of  a  body  of  revolution 


One  of  the  quantities  that  is  studied  is  the  product 


(=i  -  t) 


K,  where 


s-i - is  the  distance  of  the  nose  of  the  shock  wave  from  the  nose  of 

x  2 

the  sphere  and  K  is  the  curvature  of  the  shock  wave  at  its  nose.  Two 


approximations  to  the  variation  of  |^S]_  -  — jK  with  Mach  number  are 

derived.  See  equations  (29)  and  (30)  of  reference  5*  These  equations 
have  been  used  to  plot  the  curves  shown  in  figure  9-  These  theoretical 
results  can  be  compared  with  the  present  experimental  results.  The 
curvature  K  is  easily  shown  from  equation  (2)  to  be  Experi¬ 

mental  values  of  the  product  /sj_  -  — \k  have  been  computed  from  the 

\  27 

values  of  c2  in  table  II  and  the  values  of  s^/d  in  figure  6.  These 
experimental  values  are  shown  in  figure  9 • 


Location  of  sonic  line.-  For  calculation  of  the  pressure  drag  by 
the  method  of  reference  1,  a  knowledge  of  at  least  the  approximate 
location  of  the  sonic  line,  as  well  as  the  shock  wave,  is  needed.  The 
two  ends  of  the  sonic  line  can  be  located  fairly  accurately  from  the 
results  given  in  the  present  paper.  The  location  of  the  end  of  the 
sonic  line  that  lies  on  the  shock  wave  can  be  found  from  the  equation 
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of  the  shock  wave,  equation  (2),  and  the  simple  theory  of  oblique  shock 
waves,  which  gives  the  angle  of  the  shock  wave  for  which  the  Mach  number 
immediately  behind  the  shock  wave  is  unity. 

The  position  of  the  other  end  of  the  sonic  line,  at  the  surface 
of  the  sphere,  can  be  obtained  approximately  from  information  that  has 
been  obtained  in  connection  with  another  investigation.  Interferograms 
of  the  flow  around  spheres  at  Mach  numbers  of  1.30  and  1.62  were  analyzed 
sufficiently  to  locate  the  end  of  the  sonic  line.  The  sonic  line  was 
found  to  intersect  the  sphere  48°  from  the  nose  of  the  sphere  at  a  Mach 
number  of  1.30  and  52°  from  the  nose  at  a  Mach  number  of  1.62.  Over  at 
least  most  of  the  range  of  Mach  number  covered  in  the  present  paper, 
therefore,  the  end  of  the  sonic  line  is  close  to  50°  from  the  nose  of 
the  sphere. 


Cones 

A  set  of  conical-tip  cylinders  was  fabricated  from  steel  rods  of 
various  diameters  by  machining  one  end  to  the  shape  of  cones  with 
various  vertex  angles.  The  circle  of  intersection  of  the  conical  part 
with  the  cylindrical  part  is  referred  to  as  the  shoulder  of  the  cone. 
Interferograms  of  the  flow  around  these  cones  at  various  Mach  numbers 
were  taken  and  are  shown  as  figures  10  to  15. 

Axial  distance  of  shock  waves.-  The  axial  distance  s2  of  the 
shock  wave  from  the  plane  of  the  shoulder  was  measured  on  the  interferograms. 
The  ratio  s2/d  is  shown  in  figure  l6  as  a  function  of  Mach  number.  A 
simple  equation  that  describes  the  variation  of  Sg/d  with  Mach  number 
over  the  range  of  the  experiments  was  obtained  from  the  data.  The 
equation  is 


—  =  0.48(m  -  l)"1/2  (5) 

d 

and  was  used  to  plot  the  curve  shown  in  figure  l6. 

Although  a  wide  range  of  vertex  angle  was  used,  S2/d  does  not 

vary  much  with  vertex  angle  at  a  given  Mach  number.  The  variation 
in  s2/d  with  angle  does  not  exceed  about  10  percent.  This  result 
was  true  even  when  the  vertex  semiangle  was  increased  to  90°  and  the 
end  of  the  cylinder  thus  became  flat.  This  result  agrees  with  Busemann1  s 
result  (reference  8),  which  is  stated  for  two-dimensional  flow  but  which 
should  be  true  also  for  three-dimensional  flow,  that  the  shoulder  is  the 
most  important  part  of  a  body  in  determining  the  detached  shock  wave 
and  that  the  nose  is  unimportant  in  the  determination  of  the  location 
of  the  detached  shock  wave. 
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The  result  that  the  cone  angle  had  little  effect  on  the  value  of 
S2/d  was  not  true  when  the  cone  angle  was  decreased  enough  for  the  shock 
wave  to  approach  attachment  to  the  tip  of  the  cone,  A  typical  case  when 
the  shock  wave  is  nearly  attached  is  shown  in  figure  15.  The  distance 
S2  given  by  equation  (5)  does  not  apply  in  this  case. 

Equations  of  shock  waves.-  Measured  values  of  the  coordinates  of  the 
shock  waves  are  given  in  table  III.  In  a  manner  similar  to  that  used 
for  the  spheres,  a  power  series 

X  =  c2Y2  +  c^Y4 

was  fitted  to  the  shock  waves  on  the  cones.  Computed  values  of  the 
coefficients  C2  and  c^  and  of  the  root-mean-square  deviation  5 
are  listed  in  table  IV. 

Location  of  sonic  line.-  As  was  the  case  for  spheres,  the  two  ends 
of  the  sonic  line  can  be  located.  The  end  that  lies  on  the  shock  wave 
can  be  found  from  the  theory  of  oblique  shock  waves.  The  other  end 
lies  at  the  shoulder  of  the  body  and  is  perpendicular  to  the  surface  of 
the  conical  tip  (reference  2). 


CONCLUDING  REMARKS 


Data  have  been  obtained  on  the  shape  and  the  location  of  detached 
shock  waves  on  spheres  and  cones.  These  data  have  been  correlated  to 
give  equations  that  describe  the  shock  waves.  Information  about  the 
location  of  the  sonic  line  has  also  been  given.  With  the  boundaries 
of  the  subsonic  region  behind  the  detached  shock  wave  thus  known,  the 
pressure  distribution  and  the  pressure  drag  of  the  fore  part  of  spheres 
and  cones  can  be  calculated,  as  for  example  by  the  method  of  Maccoll. 


The  experimental  data  given  herein  should  also  be  useful  for 
comparison  with  the  results  of  future  developments  in  the  theory  of 
detached  shock  waves. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 

Langley  Air  Force  Base,  Va.,  September  19,  19^9 
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TABLE  I.-  COORDINATES  OF  SHOCK  WAVES  ON  SPHERES 
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Figure  1.-  Interferogram  of  flow  ar 
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Figure  2. - 


L-62131 

Interferogram  of  flow  around  sphere  at  Mach  number  I.30. 
d  =  i  inch. 
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Figure  5*-  Interferograms  of  flow  around  spheres  at  Mach  number  1.8l 
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Figure  8.-  Variation  of  coefficients  c2  and  c^  with  M  for  spheres. 
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Figure  9*“  Comparison  of  experimental  and  theoretical  results  on  product  of  shock-wave  distance 

curvature . 
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Figure  11. -  Interferograms  of  flow  around  cones  at  Mach  number  1.30. 
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(b)  6  =  45°;  d  =  -2^  inch. 
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Figure  11.-  Continued.  M  =  1.30. 
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(a)  e  =  35°;  d  =  ^  inch. 
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Figure  12.-  Interferograms  of  flow  around  cones  at  Mach  number  1.37 
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Figure  12.-  Continued.  M  =  1.37* 
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Figure  13 . -  Continued.  M  =  1.62. 
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Figure  l4.-  Continued 
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Figure  l6.-  Variation  of 
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